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SPECTRAL  ANALYSIS  OF  THE  ELATINA  VARVE  SERIES 

R.N.  Bracewell 

Space,  Telecommunications  and  Radioscience  Laboratory,  Stanford,  California  94 305-4035,  U.S.A. 


Abstract.  The  Elatina  formation  in  South  Australia,  which  provides  a  rich  fossil  recuid 
of  presumptive  solar  activity  in  the  late  Precambrian,  is  of  great  potential  significance  for  the 
physics  of  the  sun  because  it  contains  laminae  grouped  in  cycles  of  about  12,  an  appearance 
suggestive  of  the  solar  cycle.  The  actual  spectrum  of  the  lamina-thickness  series  is  rather 
complex,  20  or  more  spectral  lines  having  been  recognized  by  Fourier  analysis.  It  is  shown 
how  these  numerous  lines  arise  as  combination  frequencies,  from  a  much  simpler  intrinsic 
spectrum,  by  rectification.  Optical  studies  of  the  sun  have  shown  that  there  is  a  magnetic 
polarity  reversal  on  the  sun  every  11  years  approximately,  but  terrestrial  consequences  of 
solar  activity,  for  example  in  the  ozonosphere  or  ionosphere,  do  not  respond  to  solar  mag¬ 
netism;  thus  the  negative-going  semicycles  of  the  full  magnetic  cycle  are  in  effect  rectified. 
Application  of  this  knowledge  to  the  Elatina  formation  shows  that  derectification  simplifies 
the  spectrum  of  the  lamina-thickness  series  in  exactly  the  way  that  one  would  expect  if  the 
solar  cycle  were  at  w'ork  here  also.  Zig-zag  effect,  an  alternation  of  cycle  thickness,  is  taken 
to  be  due,  not  to  a  beat  phenomenon,  but  to  rectification  in  the  presence  of  a  weak  345-year 
oscillation;  subtraction  of  this  oscillation  after  dei edification  is  essential  to  the  simplifying 
procedure.  The  fundamental  period  is  established  at  a  new  sharper  value  of  23.7  ±  0.2  years 
as  compfured  with  the  looser  22.2  ±  l.S  years  for  the  modern  sunspot  series.  This  paper 
treats  the  laminae  as  varves  laid  down  yearly  and  modulated  in  thickness  in  accordance 
with  the  late  Precambrian  sunspot  activity  for  the  year  of  deposition.  Since  the  difference 
between  23.7  and  22.2  is  less  than  a  standard  deviation  it  is  premature  to  speculate  that 
the  sunspot  cycle  period  has  undergone  secular  change:  indeed  the  possibility  that  the  solar 
oscillator  has  been  secularly  stable  is  not  ruled  out.  The  high  Q  now  demonstrated  for  the 
varve  oscillator  (around  120  compared  with  a  previous  value  of  12)  weakens  support  tor 
that  part  of  solar  dynamo  theory  that  ascribes  the  solar  cycle  to  a  self-sustaining  relaxation 
oscillation;  conversely,  the  evidence  fen-  an  internal  solar  clock  mechanism  is  strengthened. 
A  wave  propagation  zone  intervening  between  the  clock  and  the  solar  surface  could  produce 
the  intrinsic  spectrum. 


1.  Introduction 

Sedimentary  laminae  of  light-coloured  silt  or  fine  sand  in  a  dark  clay  matrix,  found  in  the 
6S0  million  year  old  Elatina  formation  at  Pichi  Richi  pass,  South  Australia,  contain  a  record 
that  has  been  interpreted  as  annual  solar  activity  in  late  Precambrian  time  (Williams  1981, 
19S3,  1985,  1986,  1987;  Williams  and  Tonkins  1985).  This  paper  is  one  of  a  series  (Bracewell, 
1985,  1986c,  1988a, b,c,  Bracewell  and  Williams  1986)  stimulated  by  the  varve  interpretation. 
Figure  1  presents  the  varve-thickness  series  for  1337  years  that  has  been  studied  in  previous 
papers  (Williams,  19S5;  Williams  and  Sonett,  1985;  Bracewell  and  Williams,  19S6;  Sonett 
and  Williams,  1987;  Bracewell,  19S8a).  In  this  informative  but  complicated  graph  of  the 
actual  data  one  sees  principal  peaks  P\  of  the  Elatina  cycle  spaced  314  years,  subsidiary 
peaks  A,  and  episodes  of  antiphased  zig-zag  effect  (alternating  peak  value)  at  A  and  B. 
Zig-zag  effect  presents  itself  to  the  eye  as  an  alternation  in  thickness  of  the  groups  of  about 
12  laminae  comprosing  the  cycles  that  are  identified  with  the  solar  cycle  under  the  varve 
interpretation.  An  alternative  interpretation  offered  by  Zahnle  and  Walker  (19S7)  ascribes 
zig-zag  effect  to  tidal  influence  but  cannot  apply  to  the  Elatina  formation  because  the  two- 
frequency  beat  mechanism  proposed  calls  for  pronounced  modulation  of  the  thinnest  laminae 
of  each  cycle. 

Following  a  preliminary  study  (Williams,  1985),  more  detailed  spectral  analysis  of  the 
1337  varves  was  carried  out  by  Sonett  and  Williams  (1987),  who,  in  a  major  investigation, 
not  only  obtained  the  discrete  Fourier  transform  but  also  a  maximum  entropy  estimate  of 
the  spectrum,  a  least-squares  periodogram,  and  a  comparison  between  the  two  halves  of 
the  sequence.  Table  I  summarizes  the  frequencies,  periods  and  amplitudes  of  the  40  spectral 
peaks  recorded  in  their  Tables  I  la  and  lib.  For  case  of  comparison  with  the  present  paper,  the 
amplitudes  presented  are  the  square  roots  of  the  original  power  spectral  densities,  normalized 
to  make  the  largest  value  equal  to  100.  Fig.  2(a)  shows  the  spectral  lines  from  Table  1,  those 
from  the  upper  half  of  the  varve  sequence  being  drawn  upward  and  those  for  the  lower  half 
downward. 

Clearly  the  actual  spectrum  is  a  complex  one  that  is  not  easy  to  absorb.  The  purpose 
of  this  paper  is  to  present  a  simpler  structure,  or  intrinsic  spectrum  [Figure  2(b)]  that  will 
be  uncovered  by  appropriate  data  analysis. 


•i 


TABLE  I 


Spectral  peaks  found  in  the  upper  and  lower  halves  of  the  varve-thickness  series  by  Sonett 


and  Williams  (19S7). 


Upper  half 

Lower  half 

Frequency  Period 

c/a  yr 

Amplitude 

Frequency  Period 

c/a  yr 

Amplitude 

0.003 

334 

73 

0.003 

332 

87 

0.006 

167 

67 

0.006 

165 

51 

0.010 

99.2 

19 

0.011 

SS.3 

32 

0.016 

62 

S 

0.036 

27.9 

22 

0.026 

38.7 

6 

0.039 

25.6 

3S 

0.033 

30.7 

19 

0.042 

23.9 

35 

0.039 

25.7 

51 

0.045 

22.3 

40 

0.041 

24.3 

17 

0.048 

21 

16 

0.045 

22.3 

59 

0.051 

19.6 

17 

0.048 

20.9 

13 

0.055 

18.0 

17 

0.051 

19.6 

31 

0.058 

17.2 

16 

0.062 

16.2 

8 

0.061 

16.3 

22 

0.065 

15.4 

14 

0.065 

15.4 

15 

0.073 

13.8 

35 

0.068 

14.7 

29 

0.077 

13.0 

67 

0.070 

14.2 

41 

0.081 

12.4 

56 

0.075 

13.4 

51 

0.084 

11.9 

100 

0.081 

12.4 

98 

0.087 

11.5 

42 

0.084 

11.9 

SO 

0.090 

11.1 

SS 

0.090 

11.1 

79 

0.097 

10.3 

46 

0.096 

10.4 

52 
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2.  The  Centralized  Series  C(t)  and  the  Modulation  Envelope 

Among  phenomena  that  are  apparent  from  Figure  1  is  a  pronounced  modulation  of  the 
layer  thickness,  referred  to  as  the  Elatina  cycle,  that  has  a  period  Te  of  about  314  years.  The 
modulation  envelope  was  extracted  by  Williams  (19S3,  1985),  and  redetermined  by  Bracewell 
and  Williams  (1986)  as  follows.  The  varve  series  that  gives  the  thickness  in  pm  of  the 
varve  for  year  t  was  first  derectified  to  yield  the  alternating  series  V±(t)  defined  by 

V'(f)  =  Vb+  1 14(0  1  • 

In  this  equation  Vo  is  a  numerically  determinate  threshold  of  70  pm  interpretable  as  the 
thickness  of  silt  that  would  be  deposited  in  a  year  when  the  sunspot  number  was  zero. 
After  Vo  was  subtracted  from  V(t),  the  sign  of  the  difference  was  reversed  at  the  times  t} 
of  minima  and  the  value  V± (t})  was  put  to  zero.  In  appropriate  crises  the  sign  reversal  was 
taken  between  years. 

Interpretation  of  zig-zag  effect.  A  key  feature  of  the  present  procedure  is  the  removal 
of  zig-zag  or  sawtooth  effect  (Williams,  19S5).  It  had  been  mentioned  (Williams,  19S5) 
that  the  sawtooth  pattern  (a  line-segment  graph  of  cycle  thickness  versus  cycle  number)  is 
mathematically  describable  as  a  beat  between  two  oscillations  with  periods  near  24  years. 
Sonett  and  Williams  (I9S5)  suggested  that  a  beat  between  periods  22.4  and  25.8  years 
might  be  responsible,  but  examination  of  Figure  2(a)  warns  that  the  two  periods  mentioned 
be  taken  in  context.  When  all  the  components  of  the  full  24-year  multiplet  beat  together 
the  result  is  a  waveform  not  possessing  the  necks  that  Williams  (19S5)  described. 

The  interpretation  (Bracewell,  1985)  of  sawtooth  effeci  to  be  adopted  here  is  illustrated 
in  Figure  3.  The  upper  part  shows  an  oscillation  of  short  period  T_\i  (representing  the 
magnetic  cycle)  superimposed  on  a  weaker  oscillation  of  (longer)  period  T\j  =  345  years; 
when  the  composite  waveform  is  rectified  (lower  part)  a  double-branched  upper  envelope 
is  produced  within  which  zig-zag  effect  is  seen.  The  double  upper  envelope  crosses  itself 
(or  passes  through  necks  Ar)  at  intervals  7y/2,  and  at  each  such  neck  the  phase  of  the 
zig-zag  effect  reverses  by  ISO0.  This  elementary  explanation  of  zig-zag  effect  shows  that 
the  upper  envelope  of  the  example  has  two  sinusoidal  branches  that  intersect,  and  that  the 
rectified  waveform,  representing  varve  thickness,  rises  up  from  a  single-valued  level  base  to 
the  double- valued  upper  envelope,  rising  alternately  to  each  branch. 

A  t'w-'dement  beat  pattern  [such  as  a)  sin(2~f/22.4)+U2  sin(2~f/25.S)]  docs  not  rise  from 
a  defined  base  level  but  shuttles  quasisinusoidally  back  and  forth  between  a  double  upper 
envelope  and  a  lower  envelope  that  is  not  level  and  is  also  double- valued.  The  envelopes  are 
not  sinusoidal  and  do  not  intersect,  unless  a\  =  a 2. 


When  one  sums  the  full  multiple  Ea,  sin(2-f/P,)  comprising  approximately  seven  com¬ 
ponents  with  amplitude  a ,  and  period  P,  around  24  years,  the  resultant  sum  describes  the 
Elatina  cycle,  as  to  amplitude,  and  through  modulation  of  the  times  of  zero  crossing  expresses 
phase  and  frequency  modulation  associated  with  the  Elatina  cycle.  Thus,  beats  between  two 
adjacent  frequencies  are  not  the  source  of  the  zig-zag  effect. 

An  explanation  of  zig-zag  effect  in  terms  of  a  beat  between  two  sinusoids  given  by  Zahnle 
and  Walker  (19S7)  in  support  of  a  lunar  theory  of  the  Elatina  varves  also  fails  to  reproduce 
the  single- valued  level  lower  envelope  and  therefore  cannot  be  correct.  The  present  paper 
takes  zig-zag  effect  as  evidence  for  the  existence  of  a  weak  long  period  oscillation  U{t)  of 
period  Ty,  to  be  referred  to  as  the  Elatina  undulation. 

Derivation  of  the  centralized  series.  After  derectification,  the  Elatina  undulation  £/(/), 
which  has  a  period  Ty  of  approximately  345  years,  was  determined  with  the  use  of  a  specific 
low-pass  filter  and  subtracted  from  V± (<)  to  establish  a  centralized  undulation-free  oscillation 
C(i).  Thus 

C(t)=v±(t)-U(t). 

The  preparation  of  the  data  makes  C(t).  shown  in  Figure  4.  easier  to  assimilate  than  the 
graph  of  V(t)  in  Figure  1.  As  a  further  aid  to  comprehension,  the  starting  points  of  the  four 
segments  shown  advance  314  years  at  a  time.  Not  only  is  the  modulation  envelope  much 
more  apparent  to  the  eye,  but  so  is  the  degree  of  repeatability.  Adjacent  semicycles  with 
contrasting  shapes  are  conspicuous;  the  sharp  negative  semicycles  at  years  62  and  S5  result 
from  a  near  coincidence  of  a  maximum  of  the  Elatina  cycle  with  a  negative  maximum  of 
the  undulation.  The  nonlinearity  described  in  earlier  papers  (e.g.  Bracewell,  19SSb)  acts 
to  produce  successive  sharp  and  blunt  semicycles.  Variations  in  semicycle  duration  are  also 
easily  seen. 

Elatina  cycle  envelope.  Summing  the  absolute  values  of  C{t )  over  a  double  cycle  (the 
varve  equivalent  of  the  solar  magnetic  cycle),  i.e.  from  each  t}  to  t}+ 2,  gave  a  modulation 
envelope  ( Figure  5)  that  essentially  confirmed  the  earlier  description  given  by  Williams  (1983, 
1955).  This  envelope  method,  apart  from  the  data  conditioning,  is  related  to  an  envelope 
method  described  by  Dicke  (1978). 

The  shape  of  the  modulation  envelope  is  not  simple;  it  exhibits  a  principal  peak  Pi  at 
year  GO,  a  secondary  peak  P2  at  year  27S  of  relative  strength  about  0.6,  a  broad,  principal 
minimum  Ml  at  year  1S1  that  descends  to  about  0.4  of  the  peak  value  at  Pi,  and  a  narrow 
secondary  minimum  M2  at  year  3  that  descends  to  about  0.5  of  PI.  In  addition  there  is  a 
weak  but  noticeable  wiggle  W  between  PI  and  Ml  at  about  year  102. 


Fourier  analysis  of  the  averaged  envelope  f?£(t )  gives  the  results  shown  in  Table  II  for 
the  coefficients  a„,  bn ,  cn  defined  by 

„  ,  v-'  f  .2 Trnts  ,  .  ,2irnt' 

E£{t)  =  1  +  2_^  [an  cos(-^— )  +  bn 

and  cn  =  +  b^.  The  time  origin  is  chosen  so  that  the  principal  peak  P\  is  at  year  60 

in  Figure  5,  corresponding  to  regularly  spaced  varve  years  60,  374,  6S8,  1002  and  1316  in 
Figure  1. 

TABLE  II 


Fourier  coefficients  for  the  Elatina  modulation  envelope  £?£(i). 


n 

1 

o 

3 

4 

5 

6 

7 

0.133 

-0.190 

-0.099 

-0.034 

0.022 

0.008 

0.002 

bn 

0.331 

-0.007 

-0.052 

-0.064 

-0.019 

0.009 

0.008 

0.356 

0.190 

0.112 

0.072 

0.029 

0.013 

0.008 

3.  Spectral  Analysis 

The  centralized  undulation-free  series  C(t)  was  extended  with  trailing  zeros  to  a  length 
of  2048  elements  and  Fourier  analyzed  by  means  of  the  Hartley  transform  (Bracewell,  19S6b). 
The  result  is  shown  in  Figure  6.  There  is  a  fundamental  band  around  a  24-year  period  and 
a  noticeable  8-year  third  harmonic  band.  There  is  no  zero- frequency  component  and  little 
at  the  second  harmonic  (period  12  years).  However,  the  fundamental  contains  a  series  of 
sidebands  reflecting  the  known  spectrum  of  the  modulation  envelope. 

If  one  added  the  345-year  undulation  back  into  C(t),  and  rectified,  and  then  performed 
a  Fourier  analysis  of  |  C(t)  +  U(t)  |  a  very  much  more  complicated  spectrum  would  result 
that  included  a  long  series  of  harmonics  of  a  12-year  fundamental  with  sidebands,  harmonics 
of  a  172. 5-year  fundamental,  a  band  extending  from  zero  frequency,  and  numerous  sum  and 
difference  frequencies.  This  complicated  spectrum  is  essentially  the  one  determined  by  Sonett 
and  Williams  (1987).  Many  of  the  lines  of  the  complicated  spectrum  are  simply  artefacts  of 
rectification  and  are  without  independent  importance;  conversely,  the  third  harmonic,  which 
is  important,  is  cancelled  out  by  rectification,  and  does  not  appear  in  the  Sonett  and  Williams 
spectrum  (which  will  be  referred  to  as  the  actual  spectrum),  because  its  sign  is  reversed  at  1.5 
cycle  intervals.  As  a  matter  of  history,  the  third  harmonic  was  discovered  in  connection  with 
sunspots  by  Cole  (1973),  who  was  the  first  to  apply  the  concept  of  derectification  (Bracewell 
1953)  to  Fourier  analysis  of  the  sunspot  number  series. 

Clearly  the  spectrum  of  Figure  6  describes  the  data  rather  economically.  A  detailed  dis¬ 
cussion  is  given  in  Section  6;  meanwhile  the  principal  spectral  peaks  of  C{t)  in  the  frequency 
range  from  zero  to  0.1  c/a  are  shown  in  Table  III  and  have  been  graphed  in  Figure  2(b). 
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The  iine  at  345  yr  is  not  part  of  the  spectrum  of  C(t),  but  has  been  indicated  for  reference. 
Application  of  the  idea  of  derectification  shows  that  all  the  previously  reported  spectral  lines 
derive  from  a  much  simpler  intrinsic  spectrum  comprising  (a)  a  23.7-year  carrier  accompa¬ 
nied  by  a  few  sidebands  associated  with  the  314-year  Elatina  amplitude  modulation  and  (b) 
a  weak  line  at  345  years.  That  all  the  detail  should  be  explicable  in  terms  of  derectifica¬ 
tion  was  inevitable;  what  was  not  inevitable  is  that  the  new  description  should  prove  to  be 
notably  simple,  a  result  that  confirms  the  physical  significance  of  this  unravelment. 

TABLE  III 

Spectral  peaks  of  C(t). 


Frequency 

c/a 

Period 

yr 

Amplitude 

0.029 

34.1 

11 

0.033 

30.6 

14 

0.036 

28.1 

16 

0.039 

25.6 

56 

0.042 

23.8 

100 

0.045 

22.0 

20 

0.048 

20.7 

46 

0.052 

19.3 

10 

0.055 

1ST 

10 

If  the  spectrum  of  C(t)  in  Figure  6  is  compared  with  Figure  2(a),  it  will  be  seen  that 
the  content  of  the  frequency  band  from  zero  to  0.02  c/a  has  disappeared  as  a  direct  result  of 
derectification.  There  is  nothing  in  C(t)  at  345  years,  or  submultiples,  because  the  undulation 
was  removed  in  advance.  The  weak  band  occurring  in  the  actual  spectrum  around  the  24- 
year  period  has  become  the  principal  feature  in  the  whole  spectrum  of  C(t),  while  the  12-year 
band  in  the  actual  spectrum  has  become  insignificant.  The  period  of  the  principal  line  is 
23.8  years,  but  this  coarse  value,  imposed  by  the  discrete  frequency  step  of  1/20  IS  c/a,  will 
now  be  refined. 

4.  The  Fundamental  Period 

Hie  fundamental  period  of  the  varve  series  has  been  previously  taken  as  approximately 
12  years.  Just  as  the  11-year  sunspot  cycle  series  derives  from  a  22-year  magnetic  cycle,  so 
also  the  varves  could  be  treated  as  deriving  from  a  24  year  cycle.  We  understand  clearly  that 
the  solar  cycle  variations  in  solar  ultraviolet  radiation,  in  terrestrial  ozone  concentration,  and 
in  ionospheric  electron  density,  do  not  reverse  when  the  magnetism  reverses.  The  magnetic 


cycle  has  been  rectified  before  these  effects  of  solar  activity  are  produced,  and  the  same  would 
be  true  of  any  variations  in  varve  thickness  attributable  to  the  solar  cycle.  If  we  allow  for  the 
frequency  doubling  due  to  rectification,  the  12-year  period  would  derive  from  a  double-cycle 
period  in  the  neighborhood  of  24  years.  A  precise  value  would  be  important  in  case  there  is 
a  difference  between  the  Precambrian  varve  period  and  the  modern  11-year  sunspot  period. 
A  significant  difference  could  have  a  bearing  on  the  theory  of  the  solar  activity  cycle,  unless 
of  course  the  solstitial  year  has  changed  in  700  million  years. 

Theoretical  reasons  have  been  adduced  for  thinking  that  the  sunspot  period  should  be 
lengthening,  not  shortening  (New  Scientist,  19S7;  Noyes  et  al.,  19S4)  and  used  against  the 
solar  interpretation  of  the  varves.  But  the  deduction  is  from  dynamo  theory  which  has  not 
in  fact  been  established  as  the  cause  of  the  sunspot  cycle  (Gilman,  19S1).  On  the  contrary, 
the  high  Q  found  below  presents  a  new  difficulty  for  relaxation  mechanisms  such  as  proposed 
by  the  dynamo  theory  of  the  solar  cycle. 

A  nominal  value  of  24.0  years  for  the  carrier  period  of  C{t )  results  from  noting  that  the 
first  downcross  is  at  year  6  and  the  last  at  year  1327,  which  is  55  downcrosses  later.  The 
calculation  is  sensitive  to  the  two  critical  zero  crossings  and  does  not  take  the  record  as  a 
whole  into  account.  The  value  24.0  ±  1.75  obtained  by  Eracewell  and  Williams  (19SG)  from 
a  histogram  of  cycle  durations  defined  by  zero  crossings  is,  in  effect,  defined  by  just  the  first 
and  last  downcrosses. 

A  different  approach  would  be  to  locate  the  centroid  of  the  part  of  the  spectrum  in  Figure 
2(b)  that  is  concentrated  in  the  vicinity  of  24  years.  As  can  be  seen,  the  complex  of  lines 
is  not  symmetrical,  possibly  because  of  the  fluctuations  that  are  characteristic  (Bracewell 
19SGa)  of  power  spectra,  but  it  is  evident  that  a  value  less  than  24  years  would  result. 

The  central  peaks  in  the  12-year  ^and  of  the  Sonett  and  Williams  spectrum  are  at  11.1, 
11.5,  11.9,  12.4,  and  13  years;  they  do  not  sharply  indicate  a  mid-frequency  because  the 
peaks  are  at  sum  frequencies,  arising  through  rectification,  of  fine  structure  in  the  24-year 
band  [Figure  2(b)],  This  fine  structure  represents  the  sidebands  that  would  be  expected  from 
the  known  spectrum  (Table  II)  of  the  Elatina  amplitude  modulation,  and  in  the  simplest 
case  would  consist  of  peaks  spaced  at  integral  multiples  of  1/314  c/a. 

The  24-year  complex  is  shown  to  an  expanded  scale  in  Figure  7.  A  string  of  short 
ticks  straddling  the  spectrum  indicates  the  expected  locations  of  the  eight  inner  sidebands 
associated  with  the  Elatina  cycle,  on  the  assumption  of  a  carrier  of  24  year  period  with 
pure  amplitude  modulation.  These  ticks  are  certainly  in  rough  agreement  with  the  spectral 
features  but  a  much  better  fit  results  when  the  carrier  period  is  adjusted  to  23.7  years. 
The  nine  vertical  arrows  rising  from  the  horizontal  axis  toward  the  spectral  peaks  indicate 
the  revised  sideband  locations;  the  arrows  agree  with  the  peaks  better  than  the  sideband 


locations  based  on  a  24-year  period.  Careful  consideration  suggests  that  it  has  been  possible 
to  achieve  a  precision  of  one  decimal,  or  a  fraction  of  the  width  of  the  central  line,  by 
involving  several  lines  of  the  spectrum  in  the  adjustment  of  the  period.  For  practical  use  it 
is  believed  that  a  value  23.7  ±  0.2  s.d  is  justified. 

The  new  precise  value  for  the  period  of  the  varve  double  cycle  is  longer  rather  than 
shorter  than  the  modern  magnetic  cycle  period,  but  the  difference  is  still  less  than  the 
standard  deviation  with  which  the  modern  period  has  been  determined. 

5.  The  Q  of  the  Solar  Oscillator 

The  above  discussion  lifts  the  Q  =  v / of  the  varve  osciilator  from  a  value  of  about 
12  (Bracewell  and  Williams,  19S6;  Bracewell,  198S)  to  about  120,  a  degree  of  sharpness  that 
works  in  favour  of  a  clock  mechanism  (Dicke,  1978:  Bracewell.  1988a. c)  for  solar  activity, 
assuming  that  the  varve  t.iickness  was  controlled  by  the  sun.  Some  variants  of  the  relaxation 
oscillation  idea  for  solar  activity  will  be  unable  to  account  for  such  a  high  Q.  It  is  true  that 
the  theory  of  t lie  interaction  between  electrically  conducting  solar  c  onvection  currents  and 
magnetic  fields  in  the  presence  of  solar  rotation  has  reached  a  Inch  pitch  of  development  over 
many  years,  but  whether  the  magnetic  fields  of  sunspots  can  be  generated  by  a  self-sustaining 
convection-rotation  dynamo,  or  whether  the  convection  layer  is  mainly  a  propagation  medium 
for  the  sunspot  fields,  may  now  hinge  on  the  ability  of  dynamo  theory  to  account  for  high 
Q  oscillation. 

When  the  12  year  period  was  first  reported  by  Williams  it  seemed  that  there  might  have 
been  a  speeding  up  of  the  solar  cycle  over  the  last  700  million  years  However,  it  was  not 
obvious  that  the  11  year  sunspot  period  was  incompatible  with  the  varve  period  because 
each  period  was  at  that  time  associated  with  a  Q  of  about  7  (Bracewell  and  Williams,  19S6). 
Thus  the  difference  in  periods  was  only  two- thirds  of  a  standard  deviation.  The  present 
comparison  is  between  22.2  ±  l.S  and  23.7  —  0.2.  The  separation  of  the  periods  is  still  less 
than  one  standard  deviation  and  it  would  be  premature  to  offer  explanations  for  a  slowing 
down  of  the  solar  oscillator.  Indeed  it  is  possible  that  further  consideration  may  tend  to  show 
that  the  solar  oscillator  has  changed  little  over  geological  time.  The  evidence  from  stellar 
activity  cycles,  while  limited,  dot's  not  show  dependence  of  stellar  cycle  period  on  rotation 
period  (Soderblom  and  Bulimias.  19SS). 

6.  Explanation  of  the  Actual  Spectrum 

Return  to  the  peaks  at  11.1.  11.5,  If  f),  12.1  and  13  years.  Let  /u  =  1/23.7  and  A  = 
1/314.  Peaks  in  the  spectrum  of  C(l)  are  expected  at  /()  -f  »A  ( 25.63.  23.70.  22.04.  20.59 
years).  Corresponding  to  the  various  sum  frequencies  2 /n -f  n  A  we  find  periods  1 1.02  (n  ~  2), 
11.42  (n  1),  11.85  ( n  =  0),  12.31  (n  =  —1),  12.82  (n  =  —2).  Comparing  these  periods  with 
Table  I  we  see  that  beating  between  the  carrier  and  stronger  sidebands  of  the  spectrum  of 
C(t)  is  evidently  sufficient  to  explain  the  peaks  in  the  12-year  hand  of  Sonett  and  Williams. 


Difference  irequencies  nA  likewise  arise  from  beating;  but  these  frequencies  are  so  low 
that  they  combine  with  the  345-year  undulation.  Consequently  one  expects  the  lowest  com¬ 
bination  frequencies  to  have  periods  350,  314,  314/2,  and  1/(1/314  +  1/350))  =  163.  This 
expectation  agrees  with  the  periods  near  333  and  166  years  in  Table  I  well  within  the  reso¬ 
lution  limit  reported. 

Consider  now  the  peaks  reported  with  periods  of  22.3  (22.3),  24.3  (23.9),  and  25.7  (25.6) 
years:  it  is  true  that  they  are  close  to  the  values  22.04.  23.7.  and  25.6  calculated  above. 
However,  rectification  of  a  sinusoid  with  period  24,  while  it  would  create  new  periods  at  12. 

S.  6 . would  eliminate  the  24-year  period.  Therefore  to  explain  the  24-vear  band  appearing 

in  the  Sonett  and  Williams  analysis  it  is  necessary  to  invoke  an  additional  phenomenon.  One 
contributing  cause  is  the  amplitude  modulation  but  more  important  is  the  undulation  which, 
in  the  presence  of  rectification,  doubles  any  frequencies  /o  +r)A  in  the  intrinsic  24-year  band 
and  creates  new  frequencies  at  /o  +  nA  -f  m/Tp,  m  ^  0.  The  low  order  combinations 
(n.m)  =  (0,  1).  (  — 1,  1),  and  (0,-1)  yield  theoretical  periods  22.2,  23.9  and  25.4,  to  be 
compared  with  22.3  (22.3),  24.3  (23.9),  and  25.7  (25.6)  reported.  These  numbers  are  in 
agreement  within  a  small  fraction  of  the  bandwidths  tabulated  by  Sonett  u.nd  Williams. 

It  appears  then,  that  all  the  complexity  of  the  actual  spectrum  established  by  Sonett 
and  Williams  and  presented  in  the  upper  part  of  Figure  2,  is  contained  within  the  relatively 
much  simpler  description  presented  in  the  lower  spectrum.  This  advance  in  understanding 
has  resulted  from  the  meticulous  analysis  of  Sonett  and  Williams  combined  with  the  idea 
of  derectification,  the  recognition  of  the  undulation  so  revealed,  and  the  subtraction  of  the 
undulation.  In  retrospect  one  could  say  that  the  complexity  seen  in  the  actual  spectrum 
results  from  applying  linear  techniques  of  analys.j  to  a  problem  embodying  essential  nonlin¬ 
earity.  This  is  a  satisfactory  finding  that  brings  us  closer  to  the  physics  of  the  problem  and 
is  in  the  spirit  of  an  earlier  attempt  at  simplification  by  Sonett  (19S2). 

It  is  noticeable  from  Figs.  6  and  7  that  there  is  a  continuum  underlying  the  clusters 
of  lines.  Tt  has  been  determined  that  this  continuum  arises  from  a  phenomenon  that  has 
been  reported  for  both  the  sunspot  and  varve  series,  namely  phase  modulation  (or  frequency 
modulation).  Both  phase  and  period  have  been  found  to  correlate  with  envelope  amplitude, 
and  frequency  modulation  of  the  varve  series  has  been  put  conspicuously  in  evidence  hr 
Sonett  and  Will  iams;  but  delicate  analysis  is  not  necessary  to  demonstrate  these  effects, 
which  are  visible  to  the  eye  on  the  open  time  scale  of  Figure  4.  Phase  modulation  broadens 
any  lines,  partially  fills  in  the  space  between  Imcs,  and  contributes  to  the  lowering  of  the 
apparent  Q.  Consequently  it  will  he  necessary  to  analyze  the  phase  modulation  in  order  to 
determine  the  presence  of  any  systematic  effect  tending  to  blur  the  apparent  stability  of  the 
oscillator  which,  int. insically,  may  prove  to  possess  an  even  higher  Q  than  appears  now. 
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7.  Third  Harmonic  in  the  Intrinsic  Spectrum 

A  broad  band  centered  on  S  years  represents  third  harmonic  of  the  frequencies  in  the 
'24-year  band.  The  third  harmonic  arises  from  a  phenomenon  that  has  been  explored  in 
detail  in  connection  with  sunspots  (Bracewell  19S8u).  When  a  sunspot  semicycle  is  strong, 
its  shape  departs  from  sinusoidal  in  the  direction  of  being  triangular;  in  other  words,  the 
profile  is  sharpened.  It  is  true  that,  in  the  case  of  a  strictly  periodic  waveform,  this  change  in 
shape  is  producible  by  adding  third  harmonic  (period  ~  8  years)  to  a  sinusoidal  fundamental. 
Consequently  it  is  not  surprising  that  third  harmonic  appears  in  the  spectrum  of  C(t ),  which 
after  all  is  quasiperiodic.  However,  the  "third  harmonic”  appears  as  a  band  around  S  years 
rather  than  as  a  sharp  line,  because  it  is  contributed  to  only  by  some  semicycles  (the  strong 
ones)  and  at  periods  ranging  from  6  to  10  years  in  accordance  with  the  variation  in  semicycle 
duration. 

The  sharp  semicycles  pointed  out  above  in  connection  with  Figure  4  show  that  the 
nonlinearity  responsible  also  affects  the  varves,  as  it  should  if  the  sun  controlled  the  varve 
thickness.  The  cause  of  the  sharpening  (Bracewell,  1985)  of  the  semicycle  shape  has  now 
been  attributed  to  an  empirically  determined  nonlinear  three-halves  law.  The  discovery  of 
this  law  led  to  the  realization  that  it  arises  from  the  sunspot  number  definition  (Bracewell 
1988b);  large  sunspot  groups,  being  emphasized  in  proportion  to  their  duration,  enhance  the 
sunspot  number  when  the  latter  is  large.  If  sunspots  were  due  to  an  underlying  physical 
influence  that  varied  strictly  sinusoidally  in  time  then,  despite  the  monochromatism,  the 
nonlinearity  built  into  the  sunspot  definition  would  give  rise  to  the  appearance  of  third 
harmonic  in  the  sunspot  number  series.  A  reasonable  approach  to  the  nonlinearity  would 
be  to  correct  for  it  before  performing  spectral  analysis.  The  spectrum  of  C(t)  would  then 
be  even  further  simplified,  consisting  mainly  of  a  24-year  fundamental  with  the  sidebands 
arising  from  the  Klatina  cycle  envelope. 

8.  Physical  Interpretation  of  the  Intrinsic  Spectrum 

The  24-year  multiplet  of  Figure  2(b),  comprising  a  carrier  and  several  sidebands,  could  in 
principle  be  attributed  entirely  to  the  oscillator.  However,  it  suffices  to  invoke  a  monochro¬ 
matic  oscillator,  because  sidebands  would  be  generated  as  the  cyclic  influence  was  trans¬ 
mitted  to  the  surface.  For  example,  if  the  transmission  mechanism  were  wave  propagation, 
frequency  modulation  would  be  imposed  by  a  moving  medium  such  as  the  convection  layer, 
circulation  times  of  the  order  of  centuries  being  indicated.  Likewise  amplitude  modulation 
would  result  from  propagation  through  the  inhomogeneous  solar  medium  because  of  inter¬ 
ference  between  direct  and  indirect  rays.  A  synthesis  of  the  spectrum  from  ray  theory  will 
be  presented  separately. 

9.  Discussion 
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In  t he  hope  of  understanding,  or  merely  predicting,  sunspot  activity,  authors  have  cre¬ 
ated  an  extensive  literature  of  numerical  analysis  of  the  sunspot  number  series,  but  progress 
hitherto  has  been  essentially  of  a  descriptive  kind.  The  discovery  of  the  varves  in  the  Elatina 
formation  at  Pichi  Richi  Pass  by  Dr.  G.E.  Williams  has  completely  altered  things.  The 
simple  fact  that  the  geological  varve  series  is  much  longer  than  the  sunspot  series  has  made 
it  possible  to  describe  the  solar  activity  oscillation  in  an  economical  way.  Much  previous 
complication  has  been  cleared  and  the  remaining  step  to  physical  interpretation  has  been  re¬ 
duced.  Artificial  nonlinearities  standing  between  the  physical  quantity  that  oscillates  and  the 
recorded  appearances  have  been  removed  to  give  a  glimpse  of  an  underlying  quasimonochro- 
matic  source. 

This  development  affects  two  main  topics  in  solar  physics:  self-sustaining  dynamo  action 
associated  with  the  convection  layer  and  solar  rotation,  and  the  theory  of  nuclear  power  gen¬ 
eration.  That  part  of  dynamo  theory  which  is  concerned  with  providing  a  convection-rotation 
source  for  the  solar  cycle  has  not  reached  consensus  (Piddington,  1975,  1978,  Gilman,  1981), 
which  becomes  understandable  if  the  source  is  below  the  convection  layer  (Lavzer  et  ah, 
1979).  Likewise,  solar  energy  theory,  even  though  mature,  may  benefit  from  modification  as 
to  the  composition  of  the  core.  If  the  sun  has  an  oscillating  inert  core,  then  neutrino  pro¬ 
duction  would  be  less  than  currently  calculated,  because  less  power  would  be  produced  by 
fusion  of  hydrogen  into  helium  and  more  by  the  CNO  cycle  of  reactions  (Bracewell,  19SSc). 
This  is  important  because  standard  theory  overestimates  neutrino  production  by  7. 9/2.0  and 
ways  of  modifying  the  theory  are  being  sought. 

There  is  reason  to  believe  that  in  a  sun  born  from  a  cosmic  mix  of  hydrogen,  helium  and 
‘race  elements,  iron  in  cosmic  abundance  would  be  soluble  in  the  hydrogenic  plasma  (Iyetomi 
and  Ichimaru,  1986).  However,  prior  sedimentation  in  the  protosolar  nebula  has  been  ron- 
iidered  (Prentice,  1983)  as  a  mechanism  for  a  hcavv-element  inert  solar  core:  heavy-element 
cometary  nuclei,  rare  today,  may  have  been  an  important  component  of  the  protosolar  neb¬ 
ula  (Prentice,  19SG;  Bailey  et  ah,  1986);  and  there  is  a  dynamical  argument  calling  for  iron 
and  silicate-derived  cores  in  the  major  planets  and  the  sun  (Cole,  1987).  If  comets,  or  come- 
•.esimals  (Donohue  et  al.,  1987),  became  an  important  constituent  of  the  interstellar  medium 
after  the  first  generation  of  supernovae  contributed  both  heavy  elements  and  ready-formed 
mains  such  as  silicon  carbide,  then  the  second  generation  stars  would  have  cores.  The  whole 
theory  of  the  internal  constitution  of  the  stars  would  then  require  revision. 

The  varve  discovery  will  also  affect  astrophysics  to  the  extent  that  the  dozens  of  stars  now 
evidencing  stellar  activity  cycles,  and  star  formation  in  general,  may  depend  for  elucidation 
on  progress  in  understanding  solar  activity. 


It  is  intriguing  to  see  that  fundamental  problems  of  astrophysics  are  being  affected  by 
the  geological  discovery  of  G.E.  Williams  and  to  consider  that  other  fossil  records  are  certain 
to  exist  and  that  they  will  contribute  to  the  development  of  a  new  resource  for  astrophysics. 
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TABLE  AND  FIGURE  CAPTIONS 


Table  I.  Spectral  peaks  found  in  the  upper  and  lower  halves  of  the  varve-thickness  series 
by  Sonett  and  Williams  (19S7). 

Table  II.  Fourier  coefficients  for  the  Elatina  modulation  envelope  Ev{t). 

Table  III.  Spectral  peaks  of  C(t). 

Figure  1.  Varve  thickness  V(t)  in  mm  versus  carve  serial  number  or  "varve  year." 

Figure  2.  (a)  The  two  varve  spectra  found  by  Sonett  and  Williams  after  dividing  the 
dat  a  V(t)  into  upper  and  lower  halves,  (b)  the  simpler  spectrum  of  the  centralized  waveform 
C(t)  prepared  by  derectification.  The  numerous  spectral  lines  above  prove  to  be  derivable 
from  the  simpler  spectrum  below. 

Figure  3.  Showing  how  zig-zag  effect,  with  a  period  Tv  =  345  years  (below)  arises 
by  rectification  of  an  alternating  quantity  sim2"//23.7)  superimposed  on  an  undulation 
f/sin(2rr</345)  that  is  relatively  weak  ( U  <  1).  The  scale  markers,  set  to  agree  with  the 
strong  peaks  at  the  first  episode  of  unequal  cycles  (A),  show  how  a  role  reversal  sets  in  for 
the  following  episode  ( B ). 

Figure  4.  The  centralized  undulation-free  oscillation  C(t)  prepared  from  V(t),  clearly 
exhibiting  the  Elatina  cycle  modulation,  whose  peaks  (P\,  P2)  and  minima  (Mi,  A/o)  are 
aligned  vertically.  The  345  year  undulation  has  been  subtracted  but  has  left  its  trace  in  the 
form  of  semicycles  sharpened  by  a  nonlinear  dependence  of  varve  thickness  on  the  solar  cycle 
oscillation.  Ordinates  are  normalized  to  a  maximum  absolute  value  of  unity  (at  year  1307). 
The  abscissa  is  varve  year  measured  from  origins  at  0,  314,  628,  and  942  in  successive  lines. 

Figure  5.  One  period  of  the  mean  amplitude-modulating  envelope  Eir(t)  that  character¬ 
izes  the  314  year  Elatina  cycle. 

Figure  6.  The  complete  amplitude  spectrum  of  the  centered  undulation-free  varve  series 
C(t).  Note  the  absence  of  components  at  and  around  zero  frequency,  the  faintness  of  the 
12-year  feature,  and  the  weak  band  around  8  years. 

Figure  7.  The  low  frequency  part  of  the  spectrum  of  C(l).  Vertical  arrows  show  the 
computed  location  of  sidebands  of  a  23.7  year  oscillation  subject  to  pure  amplitude  modu¬ 
lation  with  a  period  of  314  years.  Vertical  ticks  straddling  the  spectrum  and  referring  to  a 
24.0  year  oscillation,  do  not  fit  the  peaks  as  well. 
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